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Abstract 

We extend the phase field crystal method for nonequilibrium patterning to stochastic systems 
with external source where transient dynamics is essential. It was shown that at short time scales 
the system manifests pattern selection processes. These processes are studied by means of the 
structure function dynamics analysis. Nonequilibrium pattern-forming transitions are analyzed by 
means of numerical simulations. 
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1 Introduction 

A remarkable property of macroscopic systems is their ability to generate forms or patterns. 
Numerous efforts have been devoted to study main principles of pattern formation in last three 
decades. A considerable progress has been made through the study of problems related to fluid flow, 
solidification processes, formation of antiphase domain walls and grain boundaries, etc. The well 
known examples of pattern formation are: convective rolls in Rayleigh-Benard cells [l}j3], a Turing 
instability with spatio-temporal dynamics in chemical system s |4| , formation of patterns on gelation 
surfaces j5], noise induced patterns in excitable systems [6j[7] and formation of a semiconductor 
nanostructure [8j, noise induced and sustained patterns in reaction-diffusion systems |& 13 . Also 



pattern formation processes can be induced by an external influence for example irradiation 14 - 17 

A considerable study of microstructure transformation and patterning has been given using 
phase field theory. In such an approach a local atomic mass density field p = p(r, t) is introduced 
to describe the phases present. The standard phase field theory considers the field dynamics at 
diffusive time and length scales. Recently, a new direction in phase field theory known as phase 
field crystals method has been proposed (18 . This approach allows one to simulate materials on 
microscopic scales (defects and grain boundary formation) and effectively consider dynamics on 
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diffusive scales (solute transport). An advantage of this method is its ability to model a field in the 
solid phase that exhibits the periodic nature, in the framework of this approach different crystal- 
structure-related properties naturally arise in this model [19]. Moreover, this approach incorporates 



elastic and plastic behaviour of periodic systems 20 . Originally, the model was postulated in the 



framework of Ginzburg-Landau theory of phase transitions where a gradient expansion of the free 
energy is not limited by the first non- vanishing term. The free energy functional incorporates spatial 



derivatives of the field of higher order in a form of the Swift-Hohenberg operator 21 for spatial 
coupling that naturally describes periodic patterns due to minimization of the free energy. 

In the framework of the standard phase field crystals approach usually a diffusive atomic dy- 
namics is considered. Such slow dynamics can be observed in systems near its equilibrium states 
with an instant response for the change in the order parameter. In non-equilibrium systems it 
realizes if the time scale for the observed phenomenon is larger than the transient period (systems 
with dissipative dynamics). An approximation of the slow dynamics admits that a propagation 
speed of disturbances is infinitely large. Generally, if one considers systems out of equilibrium, 
then the transient dynamics should be taken into account to satisfy the criterion of a finite prop- 
agation speed according to physical microscopic processes. It follows, that the (fast) dynamics of 
such systems can be studied on the time scales of the order or smaller than the relaxation time. 



The modified phase field crystal method introduced in Ref. 22 includes both diffusive dynamics 
and elastic interactions, where the separation of time scales exists between diffusive and relaxation 
processes in solid. 

The study of the pattern formation within this approach was made in a deterministic (noiseless) 
limit. It was assumed that fluctuation sources (noise) can not principally change the system dynam- 
ics and its stationary states. Overdamped stochastic systems with nonconserved dynamics and the 



Swift-Hohenberg spatial coupling were well studied in last two decades (see, for example, 23-29]). 
It was shown that external fluctuations (additive or multiplicative) can induce pattern formation 
in such systems. Unfortunately, the patterning processes in stochastic systems with conserved 
dynamics where transient dynamics is essential were not discussed yet. 

In this paper we extend the mechanism of nonequilibrium patterning by consideration of stochas- 
tic systems with fast dynamics where the density filed is a conserved quantity. Starting from the 
balance equation we consider the system with ordinary thermally sustained flux related to fast 
dynamics and external influence leading to additional (athermal) atomic mixing. We assume that 
every flux has its own fluctuations. Considering transient dynamics we shall show that at early 
stages the system manifests pattern selection processes. We discuss nonequilibrium external noise- 
induced pattern- forming transitions in such systems. 

The paper is organized as follows. We introduce the model in Section [2] The short time insta- 
bility analysis and a possibility of pattern selection processes are discussed in Section [3j Analytical 
results we compare with computer simulations in Section [4] where nonequilibrium pattern-forming 
transitions are studied. Finally, we summarize our main conclusions in Section [5] 



2 Model 

Let us consider a class of extended systems described by a scalar conserved field x(r, t) obeying 
the mass conservation law J drx(r, i) = const. The field x relates to the mass density field p(r,t) 
as x = (p — po)/po, where p is a uniform reference density. A generic evolution model of the field 
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x is given by the continuity equation 



d t x = -V ■ J tot . (1) 

Here J tot is a total flux which generally depends on the time t. Let us assume that 3 t ot consists of 
a thermally sustained diffusion flux J d and an atomic mixing flux J e induced by external influence 
(leading to a structural disorder, turbulence effects, etc.), so that J tot = J^i + J e . The flux J e 
emerges as a result of atomic mixing caused for example by interactions of irradiated high-energy 
particles with atoms of the system (the ballistic flux) . 

Examining the system in real conditions we assume that every constituent of the total flux has 
both regular and stochastic components. The thermally sustained diffusion flux is assumed to be 
described by a relaxational Maxwell-Cattaneo equation |30| generalized by a stochastic contribution 
£, representing flux fluctuations. Therefore, the time evolution of the diffusion flux is governed by 
the Langevin equation of the form (30 32 



r D d t J D = -3 D - MV^ +£(r,t), (2) 
ox 

were td and M = const are the relaxation time and the atomic mobility, respectively; J- is the free 
energy functional of the system; £(r, t) is a Gaussian noise representing thermal fluctuations with 
(£(r,£)) = obeying the fluctuation dissipation relation (£(r, i)£(r', t)) — 2<j 2 M8(y — r')5(t — t'), a 2 
is the noise intensity reduced to the bath temperature T. The relaxation term reflects the memory 
effects; it is dominant at fast but finite speed of propagation vu = Mj tjj . In the case to — > the 
diffusion flux takes the "usual" form J^i ~ —MV8F/8x + £(r,i) as a result of an instant response 
with the diffusive speed vd — > oo. While the real diffusion processes occur at finite speed vd, next 
we consider the case td ^ 0. 

We study the model where the external flux J e satisfies the Fick law of the form J e = —D®Vx, 
D® is the external source induced effective diffusion coefficient. If one assumes that this external 
influence has stochastic nature (for example irradiated high-energy particles have Maxwell distribu- 
tion with spatial correlations (331) then one can put — D e + C(r, t), where D e is the regular part 
of the ballistic flux J e , and £(r, t) is its stochastic component describing fluctuations that emerge in 
collision processes in the system [ 34 . Considering the general problem for the external noise (,{Y,t) 
we adopt the following properties: 

(C(r, *)) = 0, (C(r, t)C(r', t)) = 2D e a 2 C(r - r')S(t - f ); 

ru i\ 1 ( (r-r'H (3) 

C (r - r ) = — = exp ' ' 



where d is the spatial dimension; a\ is the noise intensity; the coefficient D e standing in the 
correlator (|3| means that the stochastic part of the ballistic flux J e emerges only if D e ^ 0. Here 
r c is the correlation radius of the external fluctuations. Physically it corresponds to overlapping 
length of disturbed domains of atomic configurations emerging as a results of the external influence. 
In case r c — > we arrive at the white noise assumption with no overlapping. Considering a simplest 
case we assume that a 2 and r c are independent parameters. 

In the following study we consider the free energy in the standard form of the phase field crystal 
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T= I dr[f(x) + hx£x)\, i l l 
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where f(x) is the free energy density, C is the spatial coupling operator. Next we investigate a class 
of the systems which are able to produce periodic states due to minimization of the free energy T. 
As it was shown in Ref. [35] in order to construct the free energy functional for periodic systems one 
needs to take into account that in the lowest-order gradient expansion the coefficient of |Vx| 2 must 
be negative. The corresponding inifinite gradients in x need to be compensated by introduction the 
next-order terms, |V 2 ir| 2 . The next criterion for periodic pattern formation is a double well form 
of the free energy density f(x); we take it in the form 



/(z) = |z 2 + ^, (5) 

where e — 9 — 1 is the control parameter defined through the temperature T, counted off a critical 
mean filed value T c , 9 — T/T c . The spatial coupling operator C leading to formation of periodic 
patterns of the field x with a fixed wave-number has the form of the Swift-Hohenberg spatial 
interaction, C = + V 2 ) 2 , where go = 27r/a is the wave-number giving the minimum of J 7 , a is 
the equilibrium lattice spacing. With e < the homogeneous state x is unstable to the formation of 
a periodic structure for some values of the wave- vector q . Next, we put q = 1 for convenience. As 
it was shown in Ref. [18,35 a mathematical construction for C incorporates elastic effects^] Indeed, 
rewriting it in the form of the gradient expansion one has T — J dr \f{x) + [— j3\ \7x\ 2 + 7| V 2 a;| 2 ] ) , 
where j3 and 7 are some phenomenological constants. Next, if we put x — Asin(2-7rr/a), where 
f) = 1/tt 2 , 7 = 8/a 2 , and substitute it into J 7 , then we get J 7 /a ~ a -1 J dr f (x) - A 2 / clq + (4A 2 / afye 2 , 
where e = a — a . This free energy is minimized at a = a , where ag sets periodicity of the system. 
The term oc e 2 determines an elastic energy (the Hook law). Than, j3 and 7 relate to elastic 
constants fl8 36 37 . 



As a results the complete set of the dynamical equations describing the system under stochastic 
influence takes the form 

d tX = -VJ D + D e Ax + V(CVx), 

T D d t 3 D = -3 D -MV-+ e 
ox 

If one consider the limiting case of D e = (i.e., J e = 0), then combining equation for both the field 
x and the flux 3d one arrives at the phase field crystals model: TDdf t x + dtX — V • (MV^ + £). 
It follows that a hyperbolic transport here (the time derivative of the second order) emerges when 
relaxation processes of the diffusion flux are possible [3l]. As it was stated in previous studies this 
model allows one to capture a dynamics of the system at time and space scales related to molecular 
dynamics simulations (~ 10~ 12 s, ~ 10 _9 m) and consider dynamics at diffusion time and space 
scales. We generalize such the model by taking into account external stochastic influence with an 
assumption of an immediate response to the external disturbance at a distinct point. In such a 
case the presented model effectively takes into account microstructure transformation processes in 
an extended window for the time and space scales. 

It is known that the system with a hyperbolic transport can manifest pattern selection processes 
at early stages of the system evolution [38]. Our model has the same properties and should manifest 
pattern selection. Therefore, in further study we aim to investigate pattern selection processes in 
periodic systems described by the model (JgJ) under external influence. We shall consider the external 
noise induced ordering processes in the model where regular and random parts of the external flux 
have competing contributions into the system dynamics. 



1 An equivalent and alternative representation for the free energy functional {3| is T = J dr [f(x) + \ ([1 + V 2 ]x) 2 
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3 Analysis of pattern selection processes 



Considering the stochastic system one should note that only statistical measurable quantities of 
the stochastic filed are informative: the average (a;(r, t)) (a volume fraction of the system compo- 
nent) and the structure function Sk(i) (a Fourier transform of the two point correlation function 
(Sx(r, t)6x(r' , t)), 5x = x — (x)), where k is the wave- vector. To discuss a behaviour of both (x(r, t)) 
and Sk(t) we obtain dynamical equations for these quantities and analyze their solutions. 



3.1 Dynamics of the average 

To obtain a dynamical equation for the quantity (x) we average the system ^ over fluctuations 
and arrive at 

8 t (x) = -V(Jjj) + D e A(x) + V(CVx), 

&F\ (?) 



r D d t {3 D ) = -{J D )-VM , 

ox 

Noise correlators in the first equation can be decomposed using the Novikov theorem [39] 

The response function in r.h.s. of Eq.|8| can be computed from the formal solution of the Langevin 
equation ^ for the field x: 

^l = V(5(r-r>)Vx(r,t)). (9) 
Substituting Eq.Q into Eq.(|J, we get (34|[40] 

(CVx) = D e a 2 [C(r - r')| r=r , V 3 (x> + 
2 ( VC(r - r')| r=r ,) V 2 (x) + (V(*))V 2 C(r - r')| r=r ,] . 

It should be noted that C(r — r') takes the maximal value at r = r' that gives 

VC(r-r')| r=r ,=0; V 2 C(r - r')| r=r , < 0. (11) 

Hence, introducing notation ro(V 2 ) = e+C 2 +3x 2 with M = 1 where x^ represents the homogeneous 
state, we arrive at the system 

f d t (x) = -V(Jr>) + D e A(x) + D e a 2 (V 2 C(\r\)\ r=0 A(x) + D e a 2 C{Q)V\x) 
{T D d t (3 D ) = -(J D )-\7zu(\7 2 )(x). 

In our further study we move to the Fourier space. To tha t end let us introduce (x k (t)) = 
/ dr(a;(r, i))e lkr , (Jk(0) = / dr(J(r, t))e lkr and rewrite Eq.(12l in the form 

[ = -* k ( J ^k) - D e \k\ 2 (x k ) - D e a 2 V 2 C(r)\ r=Q k 2 (x k ) + D e a 2 C(0)k 4 (x k ) 

d /T x (13) 
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The system of ordinary differential equations (13) has an analytical solution 



Stability analysis of the homogeneous state x$ can be performed for the average (x k ) in the 



simplest way. Let us differentiate the first equation from the system ( 13 1 over the time t. Hence, 
expressing the flux Jdv. from the first equation and using time derivative of the flux from the second 
one, we finally obtain 



TD^fi- = (1 + T D D e k 2 E(k 2 )) ^-k 2 (D e E(k 2 ) + vo{k 2 )) (z k ), (14) 



where the notation ^(fc 2 ) = 1 + cr 2 (V 2 C(|r|) r= o — C(0)k 2 ) is introduced. A solution of the derived 
equation can be found in the form (a:| k |(i)) — (x^ (0)) exp(0(fc)i). Inserting this solution into 



Eq.(14l we get an expression for the phase 

m± = - 1 + T ^ efe2S(fc2) ± _L^(i + TD D e k*E(k*))* - 4r D k* (£> e 5(fc2) + ^(fc 2 )). (15) 

It is seen that the phase can have real and imaginary parts, i.e., <fi(k) — $t<fi(k) + i^s<j>(k). It follows 
that an unstable modes are possible only if 5R0(fc) + > 0. It is known that in systems with the Swift- 
Hohcnberg interaction there is the interval k cl < k < k c2 , where 3t</>(fc) + > with k c i,k c2 ^ 0. In 
other words the first unstable mode will always have finite period given by the wave-number from 
this interval. The quantity 5i</>(fc) + > has a unique peak always; its position sets the most unstable 
mode with the wave-number k m . From the expression ( |15[ ) it follows that the imaginary part of 
the phase emerges if the condition (1 + tdD e k 2 E{k 2 )) 2 < irok 2 (D e E(k 2 ) + o„(i 2 )) is satisfied. 
Therefore, the evolution of the average (x^(t)) can be characterized by decaying oscillations with 
the frequency 90(fc) and the decrement (1 + T]jD e k 2 E(k 2 ))/2T£) > 0. The domain for decaying 
perturbations is defined by the condition 1 + TDD e k 2 E(k 2 ) > 0. A domain for the stable modes is 
limited by the wave-number 



k 2 d = I 1 + <x 2 V 2 C(|r|) r=0 + ^(1 + a ^C(\r\) r = ) 2 + | . (16) 

A domain where the oscillating behaviour of (x| k |(i)) is realized can be defined by solutions fco = 
ko(0,D e ,a 2 ,r c ) of the equation 

1 - 2T D k 2 {D e E(k 2 ) + 2w v {k 2 )) + {r D D e k 2 E(k 2 )) 2 < 0. (17) 

3.2 Dynamics of the structure function 

Let us obtain the dynamical equation for the structure function S k . To that end we rewrite the 
system Q in the Fourier space: 

= -ikJ_D k - k 2 D e x k - k 2 CkX k 

d Jz5k , (18) 

td t; = -J-Dk - ik.w(k 2 )xk + £ k . 
at 

Than the equation for the structure function takes the form 

= -ik(J Dk a;_k) + ik(J_D_ k a; k ) - 2k 2 D e S k - fc 2 ((Ck£k^-k) + (C-kX-k^k))- (19) 



G 



The corresponding correlators can be obtained from the equation 

^ d(J Dk *_ k ) = _ {JduX _^ _ lkro(fc 2 )5k + { ^ x _^y m 

at 



The system of two differential equation of the first order (19 20 1 can be rewritten in the form of 
the one equation of the second order 

- fc 2 ((Ck^k^-k) + (C-k£-k£k)) - ik((£ k x_ k ) + (£-k£k)) (21) 

- k 2 T D ^-((( k X k X- k ) + (C-k^-k^k))- 

at 

Decomposing the noise correlators with the help of the Novikov theorem with 

(C-kZ_ k :r k ) = (CkX k x_k) = D e a 2 {V 2 C{\r\) r=0 - C(0)fc 2 )S k , 
(f-k^k) = (fk^-k) = «k6>, 

we arrive at the dynamical equation in the form 



(22) 



= - (1 + 2k 2 T D D e E(k 2 ))^ - 2k 2 (D e Z(k 2 ) + m(k 2 ))S k + 26k 2 
at'' at 



(23) 



where a\ = 6. It is seen that the obtained equation (23 1 admits a solution of the form 5* oc e^*^*, 
where the phase ^(fc), in general, can have real and imaginary parts, (p — ^Rip(k) + ^stp(k). In 
the spinodal decomposition theory with the hyperbolic transport (td ^= 0) |31| , where the spatial 
interaction is governed by the term |Va;| 2 the real part 3?(^(fc) + is known as an amplification rate 
R(k) = -$t<p(k)+, where S cx e"^*; the imaginary part Qip(k) is responsible for pattern selection 
processes. Next, we discuss behaviour of the phase <fi(k) in order to study main properties of 
the system dynamics, whereas the pattern selection processes we study considering the structure 
function dynamics. 

3.3 The external noise influence on the pattern selection 

Let us consider stability of the homogeneous state xq = studying real and imaginary parts 
of the phase 4>{k) (see FigjT]). From Figjlji it follows that in the simplest case D e = (solid line) 
unstable modes are limited by the interval of the wave-numbers k c \ < k < k C 2- The deterministic 
external influence (dashed line) suppresses the instability progress, resulting to shrinking the domain 
of unstable modes. However, external fluctuations act in opposite manner to the deterministic part 
of the flux J e : the domain of unstable modes extends as a 2 increases, and oscillating solutions 
emerge at large fc (see Fig(lb). Hence, we arrive at competing contributions of regular and random 
parts of the external flux. It should be noted that the wave-number responsible for the most unstable 
mode and, respectively, for the period of the formed structure depends on the temperature and 
the intensity D e . 

Let us study pattern selection processes considering dynamics of the structure function shown 
in Fig(2j It is seen that oscillations in time are observed at different values of the wave-number k. 
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Figure 1: Real (a) and imaginary (b) parts of the phase </>(fc) in the vicinity of the state Xq = 
at different values of the external flux parameters D e and a 2 . Other parameters are: 9 = 0.7, 
td = 0.5, r c = 0.5. 




Such a behaviour is absolutely predictable according to the form of the dynamical equation for the 
structure function. More interesting is the wave behavior of S versus the wave-number. The main 
peak in the dependence S(k) is responsible for the main period of the structures, whereas other 
peaks reflects selection of patterns with small periods. These peaks decay with time growth that 
means selection of the one unstable mode k ~ k m that gives the main contribution in the pattern 
formation processes. The same oscillations can be found in solutions of the equation for the average 
(x). From Figj2^ it is seen that during the system evolution the position of the main peak is shifted 
toward the most unstable mode k ~ fc m ; a width of the peak is reduced and it becomes higher. It 
means that the spatial patterns become well-defined with sharp interfaces. From the dependencies 
S(k) at a given t it follows that the external flux has a crucial influence on the patterns selection 
processes. Here the regular part of the athermal flux (D e ^ 0, a 2 = 0) suppresses the pattern 
selecting, whereas its stochastic contribution (er 2 ^ 0) results in the amplification of the magnitude 
of the structure function peaks sustaining the pattern selection processes. Let us note that the 
competition of both regular and stochastic contributions of the external flux leads to the fact that 



<S 
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Figure 3: Imaginary part of the phase <p(k) of state Xo = \J\ — 9 at tjj = 0.5, r c = 0.5, D e = 0.2, 
9 = 0.4. 

the main peak in S(k) dependence decreases at large D e and its width increases. It means that the 
corresponding patterns have more diffuse interfaces if additional (athermal) diffusion is introduced. 
However, the stochastic source of the flux J e acts in the opposite manner. The main peak of the 
structure function is shifted toward large k as the temperature 9 decreases. 

Let us consider stability of the state xq = \J\ — 9 at 9 < 1. We are interesting in pattern 
selection processes in the vicinity of the minimum of the free energy density f(x). From naive 
considerations one can say that the system should be stable in the vicinity of the free energy 
density minimum. Indeed, here the real part of the phase is negative, ^R(f>(k)± < always (not 
presented here). The imaginary part can exist at k > fco (see Fig{3|. Therefore, in the vicinity of 
xq = \/l — 9 pattern selection is possible too. Here an increase in the noise intensity a 1 results in 
a growth of the quantity fcrj. 

The structure function in the vicinity of the state xq = Vl — is shown in FigjZ] where os- 
cillations versus time variable and wave-number is clearly seen. The position of the main peak 
of the structure function (see Figj4^,) does not change in time at large t and its magnitude does 
not increase. It means that patterns are formed rapidly with sharp interfaces. From Fig|4|) one 
can see that a height of the main peak is reduced as D e increases, the wave behaviour disappears. 
From another hand stochastic contribution of the flux J e promotes an increase in the amplitude 
of the structure function oscillations and increase in the main peak height. Therefore, due to the 
stochastic external influence the spatial patterns become well-defined. 

In Fig{5]we plot the dependence of the period of patterns 2n/k max computed according to the 
wave-number k max related to the maximum of the structure function S max — S(k max ) at large 
t when the position of the main peak does not change. It is seen that the period of structures 
grows as temperature and the coefficient D e increases: thermal heating leads to formation of large 
patterns; and additional (athermal) diffusion results to the same effect due to the constant D e leads 
to renormalization of the effective temperature 9 e f = 1 — 9 + D e . At large D e period of patterns at 
the fixed 9 increases, but the stochastic external source suppresses this growth [^] 

dependencies in Fig[5] have the same form as (k m )~ 1 versus 8 and D e obtained from the condition of the 
maximum of the function 5R</»-|_(fc). 
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Figure 4: Evolution of the structure function (a) in the vicinity of the state x a — \J\ — 9 at 
a 2 = 0.2, D e — 0.1 and the structure function at different D e and a 2 at t = 2 (b); 9 — 0.4, 
taun = 1-0, r c = 0.65.) 



4 Simulations 

4.1 Discrete representation 

In order to verify the previous predictions, one can perform numerical simulations of the discrete 
model ([6| in a square two-dimensional lattice N x N of the cell size I. The partial differential 
equations ([6| in the discrete space take the form 



t~d 



dxi 
~dT 
dJi 
dt 



= -(Yfi)ij Jj + DeAijXj + (V R )ikCk(^L)klXh 

OF 



+ 6- 



(24) 



where index i labels cells, i = 1, . . . ,N 2 ; the discrete left and right operators are introduced as 
follows 



(25) 



(Vi)« - -(V*),-*, (V i ) <J -(Vfl) i / - An = p (<J<,i+i ~ 2^,, + 



For the stochastic sources the discrete correlator is of the form (£i(t)£j(t)) = 2£ 2 a 2 SijS(t — t'). In 
the limit -c 1 one has quasi-white noise with = 0, (£i(t)£,j(t')) ~ SijS(t — t'). Under these 

conditions we arrive at the stochastic process with properties given by Eq.([3|. 

Simulations are provided with the time step St = 10 -3 on the square lattice of the size 128 x 128 
with 1=1. Typical evolution of the system is shown in Fig{6]with initial conditions (x(r, t = 0)) = 0, 
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Figure 5: Dependencies of the period of patterns 2tt /k max versus temperature and coefficient D e 
and different a 2 at td — 0.5, r c — 0.5. 




Figure 6: Evolution of patterns at 9 = 0.3, D e = 0.5, a 1 — 1.0, td — 0.5, r c = 1.0. Snapshots are 
taken at t = 0, 50, 500, 1000. 



((8x) 2 ) = 0.3. It is seen that linear patterns are formed that corresponds to the standard analysis 
of the pattern formation in deterministic systems with such initial conditions |18| . 

4.2 Analysis of the wave behavior 

To prove the wave behavior of the first statistical moment and the structure function we analyze 
the dynamics of two first statistical moments. While the considered system obeys a mass conserva- 
tion law (J drx(r,t) = const, where const = 0), in our computer simulations we calculate averages 
that correspond to positive and negative values of the stochastic field x, i.e. (x) + = (N~ 2 J2i X T)> 
where xf relates to xi > 0, and = (TV -2 J^. xf), where xf relates to Xi < 0, (...) means 

average over experiments. If the ordering processes realize, then the quantities (x) + , (x)- should 
increase during the system evolution satisfying the conservation law (x) = (x) + + = 0. An- 
other criterion for an existence of an ordered state is growth of the averaged second moment 
J = (N~ 2 x i) which in pattern selection discussion is known as a convective heat; it plays a 
role of an order parameter in pattern-forming transitions theory |40j. An alternative definition is 
J(t) = ^2 k S(k,t), where S(k,t) is a spherically averaged structure function. Hence J(t) gives the 
square under the function S(k,t). Therefore, possible oscillations in J(t) should manifest the wave 
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behavior of the structure function versus time. To prove the existence of pattern selection processes 
we analyze the spherically averaged structure function calculated according to standard definition 
S(k,t) = (N k )-^ k ^ k+Ak S k (t).^ 

Results of our computer simulations indicate that the constituents of the total average (x(t))± 
and the order parameter J(t) increase manifesting oscillations (see Figj7^). Moreover, oscillations of 
(x) 4- and (x)— are realized in an antiphase manner that satisfies the conservation law. An increase 
of the order parameter J(t) means the ordering of the system; the corresponding oscillations reflects 
the time oscillations of the structure function. In Fig|7]3 we plot the spherically averaged structure 
function for nonlinear system. It is seen that the wave behavior versus k is realized. With an 
increase in the time t additional peaks disappear (not presented here). We plot the dependence at 
td = 1. At To < 1 these peaks are not well pronounced but exist. Heights of such peaks can be 
increased setting To > 1. In our case we consider a limit T£> < 1 assuming that time scales for the 
field x and the flux Jo can be commensurable. 

4.3 Noise- induced pattern- forming transitions 

As was shown above the interactions provided by the Swift-Hohenberg coupling operator lead 
to nonequilibrium pattern-forming transitions. In this subsection we aim to study these transitions 
by means of standard technique employed in the equilibrium phase transitions analysis. In the 
nonequilibrium case the relative order parameter is the steady state quantity r\ = lim^oo (J), where 

means average over large time interval. We can use an additional criterion m + = lim t _ i . 00 (x(i)) + 
to characterize the ordered state. The extensive related fluctuations are given by the quantity 
X — N~ 2 ((J 2 ) — (J) 2 )/(J) 2 playing a role of the generalized susceptibility. From these definitions 
it follows that in the disordered (homogeneous) state one has r\ = m + = 0, in the ordered state the 
order parameter takes nontrivial values (j] ^ 0). In the vicinity of the critical point (for example 
8 ~ C ) fluctuations grow and the susceptibility \ increases. 

In Fig(8^i we plot steady state order parameter rj and the average m + versus temperature and the 
parameter D e at different values of the external noise intensity a 2 . It is seen that as the temperature 
and ballistic mixing intensity D e grow the order parameter and the average m + decreases. In the 
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Figure 8: Plots of the steady state order parameter T), the average m+ and the susceptibility \ (a) at 
D e = 0.5 and phase diagram for nonequilibrium pattern-forming transitions (b). Other parameters 
are: tjj = 0.5, r c = 1.0. 



vicinity of the critical point 9 ~ 9 C the generalized susceptibility shows a well defined maximum 
at the same location where the order parameter departs from zero. It clearly seen that external 
fluctuations results in renormalization of the critical value for the control parameter e. One can 
find that the critical value for the noise-induced patterning tends to its mean filed value 6 C = 1 as 
the noise intensity increases at fixed D e that decreases 6 C . To illustrate the noise-induced shift of 
the critical point we compute the phase diagram shown in Fig|8]D. It follows that at fixed noise 
intensity a 2 the critical value 9 C decreases as the ballistic mixing intensity D e grows. One can 
see that if a 2 increases, then instability of the disordered state emerge at elevated temperatures. 
This conclusion is in agreement with the prediction made above from the liner stability analysis. 
Moreover, it should be noted that in the domain in the vicinity of 9 C where the fluctuations are 
large the spatial patterns have diffuse interfaces, whereas at out off critical temperatures 9 C spatial 
patterns are well-defined (see insertions in the phase diagram). 



5 Conclusions 

We have studied pattern selection processes in periodic stochastic systems with the hyperbolic 
transport. Considering the system with thermally sustained flux and flux of athcrmal mixing we 
discuss properties of pattern formation, selection and nonequilibrium pattern-forming transitions. 
The dynamics is studied in terms of both the first statistical moment and the structure function. 
Analytical results related to the linear stability analysis are compared with computer simulations. 
We have found that external athcrmal flux having both regular and stochastic components influ- 
ences crucially on pattern selection processes. It was shown that regular part of the external flux 
suppresses such processes, whereas it stochastic constituent promotes the pattern selection. Con- 
sidering pattern-forming transitions we have shown that external influence shifts the critical point 
of the transition where the regular and random components of the external flux act in competing 
manner. 
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